Generation of entangled tripartite states in three identical cavities 
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The generation of entanglement between three identical coupled cavities, each containing a single 
three-level atom, is studied when the cavities exchange two coherent photons and are in the N=2, 
4, and 6 manifolds, where N represents the maximum number of photons possible in any one cavity. 
The combined states of the atom and the photon in a cavity is given by a qutrit for A'^ = 2, a 
five-dimensional qudit for N = A, and a seven-dimensional qudit for A'^ = 6. The conservation of 
the operator N for the interacting three-cavity system limits the total number of tripartite states 
to only 6, 18, and 38, rather than the usual 3^ = 27, 5^ = 125, and 7^ = 343 states for N=2, 4, 
and 6, respectively. The deterministic generation of entanglement from general initially unentangled 
tripartite states is studied in the limit of large hopping strength, where all the solutions are analytic 
and given in terms of exponential functions. Several types of resulting tripartite entanglement are 
analyzed in order to obtain maximally entangled states. 

PACS numbers: 03.67.Bg, 03.67.Mn, 42.50.Pq 
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I. INTRODUCTION 



Quantum entanglement is essential for quantum infor- 
mation communication and processing protocols in quan- 
tum cryptography [l|, dense coding Q, teleportation 
3, and entanglement swapping Q, which can be used 
to realize quantum repeaters [5|. Entanglement can be 
achieved via two interacting quantum systems [61 or by 
an appropriate joint measurement of two systems [7|]. An 
important use of entanglement is the sharing of multipar- 
tite states distributed amongst several different parties 
that are separated by large distances where the perfor- 
mance of local measurements on their respective subsys- 
tems results in the global broadcasting of the outcomes 
of local measurements. The latter is used in standard 
quantum key distribution [8| and quantum secret shar- 
ing [9(1. Accordingly, generating multipartite entangle- 
ment is an important objective in experimental quan- 
tum systems. For instance, experimental realizations of 
four-photon entanglement and high-fidelity teleportation 
[lO| , entangled states of two and four trapped ions [ll[ , 
the entanglement of two Rydberg atoms in microwave 
cavity quantum electrodynamics (QED) |l3|, optically 
induced entanglement of excitons in a single quantum 
dot [131, creation of Greenberger-Horne-Zeilinger (GHZ) 
states with up to 14 qubits |14l |. and the entanglement 
of a six-photon symmetric Dicke state [15| . On the the- 
oretical side, many methods based on cavity QED have 
been proposed, for instance, a hybrid quantum correlated 
tripartite system formed by an optical cavity and a mi- 
crowave cavity [iq], generating n-qubit GHZ entangled 
states with a three-level qubit system and (n — 1) four- 
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level qubit systems in a cavity |17l|, and the generation 
of maximally entangled GHZ state as a ground state of 
a three spin system [l^. In general, the more particles 
or states can be entangled, the more clearly non-classical 
effects are exhibited and the more useful such entangled 
states are for quantum applications [19l |20j . 

In the study of three-level atoms, a transformation 
was introduced }21j whereby the three-level atom was 
reduced to a corresponding two-level atom of the Jaynes- 
Cummings type albeit with two-photon rather than 
single-photon transitions. This model has been used in 
cavity QED to generate "macroscopic" qubits Q and 
in the scattering of two coherent photons inside a one- 
dimensional coupled-resonator waveguide that operates 
as an ideal quantum switch [23| . The model has been ap- 
plied to two-photon exchange between two separate cav- 
ities [23 with each cavity containing a three-level atom 
in a cascade (or ladder) configuration |25l| and coupled 
via a two-photon hopping interaction j24| . The coupling 
for two-photon hopping would require a nonlinear media 
that niay be governed by a quantum Kerr-type inter- 
action [26|. The latter work 2^ was restricted to the 
N = 2 manifold, where TV denotes the maximum num- 
ber of photons possible in a given cavity. More recently, 
we extended our study to the dynamics of the iV = 4 
manifold and showed how the temporal development of 
the coupled two-cavity system generates maximally en- 
tangled states in both the N — 2 and A^ = 4 manifolds 
from an initially unentangled state [27j. 

In order to realize maximally entangled states, one 
must consider high Q-factor cavities and excited states 
that are long-lived. Owing to angular momentum and 
parity conservation, the simultaneous emissions of two 
photons possess increased lifetimes and thus are essen- 
tially metastable states. For instance, single-photon 
atomic transitions are of the order of 10~®s whereas the 



two-photon emissions have much longer hfetimes, of the 
order of fractions of a second. This makes the study 
of two-photon processes quite interesting as compared 
to single-photon transitions. In addition, two-photon 
spectroscopy has been recently used to probe the hy- 
bridization between a superconducting phase qubit and 
an intrinsic two-level system coupling to the qubit cir- 
cuit [28]. The detailed data on experiments allow the 
mapping out of this hybrid system, combining two co- 
herent quantum systems that are fundamentally differ- 
ent in nature [28|. The experimental realization of the 
transmon qubit, which is an improved superconducting 
charge qubit derived from the Cooper pair box, consid- 
ers two-photon transitions that could help design more 
robust quantum computers |29|. A recent study of the 
possibility of coherent reversal information transfer be- 
tween superconducting charge qubits and mesoscopic ul- 
tracold atomic ensembles coupled to a microwave copla- 
nar waveguide is based on two-photon optically excited 
Rydberg transitions-from the ground state to the Ryd- 
berg state via a nonresonant intermediate Rydberg state 
[30| . The experimental implementation of the latter pro- 
posal is currently underway to demonstrate hybridization 
of solid-state and atomic quantum devices |30| . 

In this paper, we consider a system of three identical 
cavities each with a three-level atom where the cavities 
are coupled to each other by a hopping interaction that 
exchanges two coherent photons. We consider the gen- 
eration of entangled tripartite states from initially unen- 
tangled states in the manifolds N=2, 4, and 6, which lie 
in a 6-, 18-, and 38-dimensional Hilbert space, respec- 
tively. Our system of three identical cavities consists of 
three atoms and one, two, and three pairs of photons in 
the manifolds N=2, 4, and 6 , respectively. Accordingly, 
generalized measurements |31j . by using only linear el- 
ements and particle detectors, may be possible on the 
generated entangled pure states. 



II. HAMILTONIAN 

The Hamiltonian for the three coupled cavities |27] is 
^iven by 



H = Y: [H^-^ - i^^l + Hiatal 
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where the first term on the right-hand-side represents the 
sum of the Hamiltonian for each cavity and the succeed- 
ing terms are the hopping interactions coupling the three 
cavities in a symmetric fashion. Therefore, the Hamilto- 
nian (1) is symmetric in the three identical cavities. 

The eigenvalues E^ and eigenfunctions |^i^ )'■'■' oi H^^^ 
are best given in terms of the dressed-atom representa- 
tion. One has 123 that 



!*,;)« =cos0„|e,n)«-sin0„|g,n + 2)«, (2) 

where \g}^^^ and je)''-' are the ground and the excited 
atomic states, respectively, \n)^^^ is the photon number 
eigenstate and 



cos^n = 

sin^n = 



r(n-f 2)1/2 



[n(r2 + l) + 2r2 + l]i/2 

jn + iy/^ 

[n(r2 + l) + 2r2 + l]i/2' 
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with r = gi/g2- 

The respective eigenvalues are given by 

^^ ^ / ^ En -\- Ep A 

E+^nu:{n+l)+ ' ^ '- -- 



+ i{A2 + ^h^glin + 2) + gl{n + l)]}i/2 (4) 



E~ = hLo{n + l) 



where gj {j — 1,2) are the atom-photon coupling con- 
stants in the three-level atom. Eg is the energy of the 
ground state, E,, is the energy of the excited state, and 
A = (Eg — E2) + hu = {Ee — E2) — TiLo is the detuning 
parameter of the mid level with energy £'2 of the three- 
level atom in a cascade (or ladder) configuration. The 
eigenenergies of the dressed states depend on n reflect- 
ing the exact treatment of intensity-dependent dynamic 
Stark shifts. 

The eigenstates l^^)*-*' are simultaneous eigenstates 
of i/(*' with eigenvalues E^, of Wq = huj{a\ai + ail — 
4g ) + [Eg + Ee)/2 with eigenvalues E' , and of iV^*) = 
ajai + (Tee — cTgg + 1 with eigenvalues n + 2. The opera- 
tor N = X]i=i A"*-'-*, associated with the three interacting 
cavities, is a constant of the motion and in this paper we 
consider the manifolds with eigenvalues A^ = 2, 4, and 
6. The set of eigenstates |*^)^*\ with n = 0, 1, 2 •••, to- 
gether with the states |g, 0)(*' = -[I'ls)^*'' and |g, 1)(') = 
|^li)'*\ where the first index refers to the ground state 
and the second to the photon-number occupation, forms 
a complete basis. 



III. N=2 MANIFOLD 

The Hilbert space of the N = 2 manifold corresponds 
to each atom in the three-cavity system being described 
by the qutrit \g,2), \g,0), and |e,0). However, owing to 
the constancy of the total operator TV = iV^^) -I- N^^^ + 
N'^^\ the space is spanned by only 6 rather than 3'^ — 27 
vectors. 

The most general uncntanglcd initial tripartite state in 
the N — 2 manifold is 



|<J+)W =sin6'„|e,n)W-f cos6'„|5,n + 2) 



(») 



1^(0)) = |g, 0)W|g, 0)(2) U, 2)(3) + 6|e, 0)(3) 



(5) 



where |ap + \b\^ = 1 and \g)^^^ and |e)*^*^ are the ground 
and the excited atomic states, respectively, and \n)'^^^ is 
the photon number eigenstate with i = 1,2,3 denoting 
the cavity. The initial state (5) is only symmetric in the 
exchange 1 -;->■ 2 and cavity 3 is in an entangled state given 
by the linear superposition of dressed states j^g )^'^^. Ow- 
ing to the symmetry of the Hamiltonian (1), the general 
solution preserves that symmetry and is given by 



m)) = \g, 0)(i)|5, 0)(2) \A{t)\g, 2) (3) + C(t)|e, O)^^) 



+ ^B{t) [\g, 0)«|.g, 2)(2) + |.9, 2) W|5, O)^^)] \g, 0)(^\ (6) 

The Schrodinger equation of the motion ih^\il}) = H\i{j) 
gives, with the aid of Eqs. (1) and (6), 

iA^ A + 2V2£,B + tan 0qC 



iB = 2V2^A + (1 + 2^)B, 



iC = tan 00 A + taiT OqC, 



(7) 



with initial conditions ^(0) = a, B{0) = 0, and C(0) = b. 
In Eq. (7), we have introduced the dimensionless time 
[{E^ — E^) cos^ Oo]t/h — > t and the dimensionless hop- 
ping coupling TiS,l [{Eq — Eq ) cos-^ 6*0] — > £, with the an- 
gle ^0 representing the mixing angle (3) for the atomic 
dressed states. The system of equations (7) gives rise to 
analytic solutions for arbitrary values of the parameters 
r and £,. However, we shall consider the case of large 
hopping strength ^ ^ 1 when the exchange of two coher- 
ent photons occurs at a much faster rate than the rate 
of atomic transitions. In that case, one has the separate 
probability conservation 



\A{t)\' + \Bit)\^^\a\^ and |C(t)p « |6p 



(8) 



Note that the amplitude C{t) associated with cavity 3 be- 
ing in the excited state is separately conserved according 
to (8). This follows directly from considering the limit of 
large hopping strength, viz, ^ ^ 1. Similarly, for the am- 
plitudes A{t) and B{t) associated with all three cavities 
being in the ground state that are separately conserved 
as indicated in (8). 

Therefore, one achieves maximal entanglement for 
times T when |^(T)p -t- |C(r)p is a minimum, that is, 
\A{T)\ is a minimum with C(0) = 6 = 0. The solution 
for the system of equations (7) is then for ^ ::!> 1 

A{t) « i (e-4'«* + 2e2'«*), 



B(i)«^(e-4'«*-e2»«*), 



Cit) = 0, 



(9) 



where the amplitudes A{t) and B{t) are periodic with 
period T — n/^. Result (9) agrees with the general re- 
sults given by (A3) in Appendix A for the corresponding 
amplitudes. Note that the state (6) becomes unentan- 
gled for times r such that ^r = 7m/3, n — 1, 2, • • •, that 
is, B{t) — 0. Our characterization of maximal entan- 
glement for the N = 2 manifold given above is quite 
consistent with the notion of geometric measure of en- 
tanglement 0-[3l|, 

Emt)))= min -\og,{mm)n, m 

|0)ePROD 

where PROD is the set of product states in the N = 2 
manifold. Now 



max 



|0)ePROD 

_ 1 
~ 9 



\{ci^m))\' = \A{t)\' + \c{tw 



(11) 



5-f4cos(6^i) 



\b? 



Therefore, the measure of entanglement < E{\il}{t))) < 
log2 9 « 3.170, where the lower bound corresponds to 
times £,T = Trl/3, / = 0, 1, 2, • • •, and the upper bound to 
times ^/Lt = (2n -I- l)7r/6 with n = 0, 1, 2, • • • and \a\ = 1. 
Note that the state (6) has an entangled two-qubit sub- 
system and, hence, does not admit Schmidt decomposi- 
tion. 

The average dwell or sojourn time the system spends 
in the entangled state governed by the amplitude 
B{t) is given by ^ J^ \B{t)\'^dt = ||a|2. Note that 
the initial state (5) possesses less symmetry than the 
Hamiltonian (1). Actually, there are only two sym- 
metric states, viz., -i^ ^p P|e,0)(i)|g,0)(2)|c,,0)(3) and 

75 Ep P\9i 2)^^^ I5, 0)(2) |g, 0)(3), where the sums are over 
even permutation, and no antisymmetric states are pos- 
sible. Both symmetric states are entangled since there 
are no unentangled, symmetric state in the N — 2 mani- 
fold. Therefore, the six-dimensional space for the N = 2 
manifold is spanned, in addition to these two symmetric 
states, by four asymmetric vectors. 



IV. N=4 MANIFOLD 

The Hilbert space of the iV = 4 manifold corresponds 
to each atom in the three-cavity system being described 
by the five-dimensional qudit \g,4:), \g,2), \g,0), |e, 0), 
and |e, 2). However, the space is spanned by only 18 vec- 
tors rather than 3^ — 125 vectors owing to the constancy 
of the operator N and so a general state is given by (Bl) 
in Appendix B. This manifold is spanned by 5 symmet- 
ric, one antisymmetric, and 12 asymmetric states, which 
treat all three cavities on the same footing. The most 
general unentangled initial states in the A^ = 4 manifold 



|V'(0)) = |g,0)«|g,0)(2) 



a|g,4)(3)-i-6|e,2) 



(3) 



(12) 



where |ap + |6p — 1, (12) is symmetric under the ex- where 
change 1 4-> 2, cavity 3 is in an entangled state, and 

m 



|^(0)) = |5,0)W 



a|ff,2)(2)+6|e,0)(2) 



c|g,2)(3)+d|e,0)(3) 

(13) ■ 

where |ap + |6p = 1, |cp + |dp = 1, and both cavities 
2 and 3 are in entangled states. For c = a and d = b, 
the initial state(13) is symmetric under the interchange 

2^3. 
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3g8«Ct _ 2e6««* + 2e-i2^?* _ 3^-4^?* 



Bi,,-^a 



gSiC* _ g6i?t _|_ g-12i«i _|_ g-4iCt 



C(„.f„ 



- e8*«* + 2e6*«* + e-i2»«* - 2e-4*«* 



(15) 



A. Initial state Eq. (12) 

The initial state (12) is symmetric under the ex- 
change of cavities 1 <-^ 2 and so the time dependent 
state {ipit)) lies in an 11-dimensional subspace rather 
than the possible 18 since A{t) = G{t), B{t) = P(i), 
D{i) = S{t), E{t) - N{t), W{t) = J(t), L{t) = C/(i), 
and M{t) = R{t) in (Bl). On the other hand, the ini- 
tial state (13) for ad ^ be ^ Q is not symmetric in the 
exchange 2 <-> 3 and so lies in the full 18-dimensional 
space. However, for ad = be the initial state (13) is 
symmetric in the exchange 2 o 3 and so the initial 
state lies also in an 11-dimensional subspace. Owing to 
the complexity of the full 18-dimcnsional vector space, 
we shall consider the limit of large hopping strength 
hS,/[{E2 —E2)cos^ 6*2] ^ 1 when the exchange of two co- 
herent photons between the cavities is much faster than 
atomic transitions. 

Consider first the solution of the Schrodinger equa- 
tions (B2)-(B7) for the initial unentangled state (12), 
that is, F(0) = a and K{Q) = 6 in (Bl). The general 
18-dimensional vector (Bl) in the A^ = 4 manifold must 
be symmetric in the interchange of cavities 1 <-> 2 and so 
the only nonzero amplitudes are given by G{t) — A{t), 
P{t) = B{t), N{t) = E{t), C{t), F{t), and K{t). There- 
fore, the solution is 



m)) « A{t) 



|5,4)W|g,0)(2) + |5,0)W|.9,4)(2) 



l5,0) 



(3) 



|ff,0)«|g,2)(2) + |5,2)W|g,0)(2) 



i3(t)|5,2)(3)+iJ(i)|e,0)(3) 



+C(i)|5,2)«|g,2)(2)|5,0)(3) 



(14) 



l.9,0)«|.g,0)(2) 



F(i)|g,4)(3)+/f(f)|e,2)(3) 



Eit) = -b 



g^iit _^ g-4iCt 



Fit) = -a 



3g8^et ^ 4g6*et ^ 2e-i2^«* + 6e-4^«* 



Kit) = -b 



2e2^«* 



-4i{i 



The amplitudes in (15) are periodic with period T = 
tt/^. Note that for any probability amplitude X(i), if 
^^^li^o = for n = 0, 1, 2, • • •, then X(i) = since 
the solutions of the dynamical equations are entire func- 
tions of t. The latter has been used in Eqs (B2)-(B7) 
to arrive at the five nonzero amplitudes given in (15) 
for solution (14). Note that 2\Eit)\^ + \Kit)\^ = \b\^ and 
2|^(f)|2 + 2|B(t)|2 + |C(i)|2 + |i^(t)|2 = |a|2 are separately 
conserved with the sum giving the overall probability of 
unity. This feature of disjoint sectors of the Hilbert space, 
determined by the number of atoms in the excited state, 
is a direct result of large hopping strength and occurs in 
all the different manifolds. 

One obtains maximal entanglement for the initial state 
(12) with a = by maximizing |_B(t)|2, or what is the 
same by minimizing [^(t)^, which occurs at a time t 
given by r^ = 7r/6, where |X(t)|2 = 1/9 and 2\EiT)\^ = 
8/9 . The measure of entanglement is logg 9 « 3.170 that 
is the same as for the N = 2 manifold albeit for different 
states. 

On the other hand, for the initial state with 6 = 
in Eq. (12), that is, all four photons are initially in 
cavity 3, there are several interesting entangled states 
obtained, for instance, by either maximizing or minimiz- 
ing |A(t)|2 + |_B(f)|2, which is the same as minimizing or 
maximizing |C(t)|2 -I- |i^(t)|2. For instance, minimization 
of |C(t)|2 + |F(i)|2 occurs at ^r ~ 0.2094,0.8378 with 
|C(t)|2 + |F(r)p Si 0.1960 resulting in a lower measure 
of entanglement log2(l/0.1960) « 2.351. The probabil- 
ities for the entangled states are 2|A(t)|2 « .2251 and 
2|B(t)|2 sa .5789. Is is interesting that at ^r = 7r/3 one 
has that Air) = B(t) = with C(7r/3) = 3V2(y3-i)/10 
and Fi-n/i) = (1 -h 3^/3i)/10. Therefore, the probability 
that all four photons remain in cavity 3 is 28% with the 
remaining 72% with cavities 1 and 2 sharing two photons 
equally. 



B. Initial state Eq. (13) 

Consider next the initial unentanglcd state (13), that 
is, L(0) = bd, M(0) = be, N{0) = ad, and P(0) = ac, 
with the time evolution of the probability amplitudes 
given by Eqs. (B2)-(B7). Now B{t) = C{t) and 
G{t) = F{t) according to Eq. (B3), D{t) = J(t) fol- 
lows from Eq. (B4), Eq.(B5) gives that E{t) = K{t), 
while T{t) = U{t) = from Eq. (B6), and Eq. (B7) 
implies that R{t) = S{t) = W{t) = 0. The state of the 
system is then, according to (Bl), given by 



-S(t)l.9,2)(i) 



|.g,0)(2)|.g,2)(3) + |5,2)(2)|.g,0)(3) 



-D{t) 



|5,2)W|e,0)(2) + |5,0)W|e,2)(2) 



l5,0) 



(3) 



+E{t) 



|ff,2)W|5,0)(^^|e,0)^^) + |g,0)W|5,0)(^^|e,2) 



i^Wl5,0)« 



|ff,0)(2)|5,4)(3) + |ff,4)(2)|g,0)(3) 



(16) 



M{t) - -be 



N{t) = -ad 



2e2^«* + g-4^«* 



2e2^«* + e-4*«* 



P{t) 
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2e«'«* + 6e6*«* + 4e-4*«* + Se-^^^^* 



21-12 and 



The solutions in (17) are periodic with period T — ir/^ 
Note that |A(t)|2 + 2|B(i)|2 + 2|i^(i)|2 + |P(i)|2 

|L(t)|2 = |6|2|d|2, 2\D{t)\' + |M(i)|2 = |6h|c , 
2|£;(i)|2 + |7V(t)|2 = |a|2|(i|2 are separately conserved with 
the sum giving the overall probability of unity. This fea- 
ture of disjoint sectors of the Hilbert space for states with 
one, two, or none of the atoms in the cavities in the ex- 
cited state is a direct result of considering the limit of 
large hopping strength. 

Consider the case with symmetry under the cavity- 
interchange 2 -f-)- 3, that is, c = a and d = 6 in Eq. (13), 
and so (17) gives that D{t) = E{t) and M{t) = N{t). 
For the initial state with a = 0, the three-cavity system 
remains in its initial state, viz., |g, 0)^^''|e, 0)^2)|g^ o)^"^-*, 
since for ^ 3> 1 the hopping interaction does not give rise 
to any temporal change. However, for the initial state 
|5,0)(i)|g,2)(2)|c,,2)(3), that is, a = c = 1 and 6 = d == 
in (13), one has 

|V(t))=A(t)|ff,4)«|g,0)(2)|g,0)(3) 



1.9,0) 



(1) 



L(t)|e, 0)(2)|e, 0)(3) + Af (t)|e, 0)(2)|5, 2)(3) 



-7V(t)|5,2)(2)|e,0)(3)+P(i)|g,2)(2)|g,2)(3) 



where 



B(t)|g,2)« |5,0)(2)|g,2)(3) + |5,2)(2)|5,0)(3) 



-F(t)|5,0)W |5,0)(2)|5,4)(3) + |5,4)(2)|5,0)(3) 



(18) 



A(t) = —ac 
^ ' 15 



- e«*«* + 2e6*«* - 2e-''*«* + e-i2»Ct 



-P(t)|g,0)«|g,2)(2)|5,2)(3), 



B{t) = -ae 



2e8*«*-3e6*«*-2e-4*«* + 3e"i2*«* 



D(t) = -be 
^ ' 3 



Eit) = -ad 



- e2*«* + e-4*«* 



- e2'«* + e-4^«* 



F{t) ^ ^ae 
^ ' 15 



gSiJt _ g6i5t ^_ g-4i5t _|_ g-12i?t 



L{t) = bd. 



(17) 



where the amplitudes are given in (17) with a = c = 1 
and 6 = d = and the normalization by |A(t)|2 + 
2|B(t)|2 + 2|F(t)|2 + |P(t)|2 ^ 1. One obtains max- 
imal entanglement for the state (18) if one minimizes 
|^(i)|2 + |P(t)|2. The same minimum occurs for ^r ~ 
.1930,0.8542,1.2402. One obtains the state where the 
respective probabilities are |A(t)|2 w 0.1070, 2|P(r)|2 k, 
0.2112, 2|P(r)|2 « 0.6060, |P(r)|2 « 0.0759. The mea- 
surement of entanglement is — log2(|A(T)|2 + |-P(t)|2) « 
2.450. 

Is is interesting that at ^r — 7r/3 one has that B{t) = 
F(t) = with A(7r/3) = V^(3 - v^i)/10 and P(7r/3) ^ 
(2 + -\/3z)/5. Therefore, the probability that the system 
is in its initial state with two photons each in cavities 
2 and 3 is 28% while the probability of four photons in 
cavity 1 is 72%. 



V. N=6 MANIFOLD 

The seven-dimensional qudit associated with the states 
of each cavity in the A^ = 6 manifold are \g,Q), |5,4), 
Ig, 2), \g,0), |e, 0), |e, 2), and |e,4), which gives rise to 
7^ = 343 states for the three-cavity system. However, 
owing to the operator N being a constant of the motion, 
the relevant Hilbert subspace for the iV = 6 manifold 
is only 38-dimensional. The dimensionality is given by 
Ix 10 -1-3x6-1-3x3-1-1x1 = 38, where the terms in 
the sum correspond to sectors a, b, c, and d, respectively 
(See Appendix C). In sector a, all three cavities are in 
the ground state and there are six photons to be shared 
pairwise amongst the three cavities. In sector b, one of 
the cavities is in the excited state, two are in the ground 
state, and there are four photons to be shared pairwise. 
In sector c, two of the cavities are in the excited state, 
only one of the cavities is in the ground state, and there 
are two photons to be shared. Finally, in sector d, all 
three cavities are in the excited state and there are no 
photons to be shared amongst the three cavities. 

The most general unentangled, initial states in the N = 
6 manifold are 



\4>m = 



a|.g,6)(i)+6|e,4)(i) 



|g,0)(2)|5,0)(3), (19) 



where |ap + |fep = 1, 



IV'(o)) = 



a|5,4)W+6|e,2)W 



c|g,2)(2)+d|e,0)(2) 



where |ap + |fep = 1, |cp + |dp = 1, and 



(20) 



IV'(O)) 



a|5,2)(i)+fe|e,0)(i) 



c|.9,2)(2)+d|e,0)(2) 



X 

(21) 



e|5,2)(^^+/|e,0) 



(3) 



where \a\^ + l&p = 1, |c|2 + jdp = l, and jep + |/|2 = 1. 

The initial state (19) is symmetric under the cavity 
interchange 2 o 3 and so |f/'(0) always remains symmet- 
rical under the interchange 2 o 3 owing to the Hamil- 
tonian (1) being symmetrical. The initial state (20) is 
strictly asymmetric. The initial state (21) may be sym- 
metric or asymmetric under the interchange of any two 
or all three cavities depending on the values of the ini- 
tial amplitudes. For instance, if a = c and b = d, then 
the solution of the Schrodinger equation will be symmet- 
ric in the interchange 1 <-> 2 for all times. If, however, 
a = c = e and b = d = f , then the solution is symmetric 
in the interchanges of all three cavities. This latter case 
is considered in Appendix D. 

In the limit of large hopping strength, ?i^ ^ [(-^'4^ ~ 
E^ ) cos^ ^4] , the vectors in the Hilbert space associ- 
ated with the sectors a, 6, c, and d are uncoupled since 
the hopping Hamiltonian in (1) only exchanges photons 



and not atomic excitations. The equations for the 38 
probability amplitudes break up into linear differential 
equations involving ten, eighteen, nine, and one ampli- 
tudes, respectively [see (C2), (C3), (C8), (C9), (Cll), 
and (C13)]. Note that for the case when all three cav- 
ities are in the excited state, the probability amplitude 
remains constant in time [see (C13)]. All these differ- 
ential equations can be solved analytically in terms of 
exponential functions for arbitrary initial states. 



A. Initial state Eq. (19) 

The initial state given by Eq. (19) lies in the non- 
overlapping sectors a and b in the large hopping limit 
Ti^ > [(-E^ — E^) cos^ Oi] presented in Appendix C. The 
time development of the system is given by a linear su- 
perposition of vectors in these two disjoint sectors. We 
consider the initial state with a = \ and so 6 = 0, that is, 
A{{)) = 1 in (CI) with the solution in the subspace of sec- 
tor a. The solution is given by (C4), where all the prob- 
ability amplitudes have analytic solutions in terms of ex- 
ponential functions and the results for A{t) and F{t) are 
given explicitly by (C5) and (C6), respectively. Maximal 
entanglement is obtained by minimizing |A(i)|^ -(- |i^(i)|^, 
which occurs for ^r « 1.7500, where |A(r)|^ -|- |i^(r)|^ « 
0.001833, 2|B(t)|2 « 0.140493, 2\E{t)\'^ « 0.055394, 
2|G(t)|2 « 0.459478, and 2\K{t)\'^ « 0.342801. This 
gives a measure of entanglement log2(546) ~ 9.1. 

On the other hand, one finds a maximum in |A(i)|^ -|- 
|i^(t)|2, which occurs for £,t « 3.0318, where \A{t)\'^ + 
|i^(T)|2 « 0.95166, |A(t)|2 « 0.89530, 2|B(t)|2 « 
0.00137, 2|£;(t)|2 « 0.02296, |F(t)|2 « 0.05637, 
2|G(t)|2 « 0.02225, and 2\K{t)\^ « 0.00176. The system 
docs not return to its initial state owing to the cigenfre- 
quencies being incommensurate. 



B. Initial state Eq. (20) 

The initial state given by Eq. (20) lies in the non- 
overlapping sectors a, 6, and c in the large hopping 
limit hS, S> [{E'l — E^)cos^ Oi] presented in Appendix 
C. Wc consider the case where the amplitudes a = 0, 
6 = c = 1, and d = in Eq. (20), that is, the 
initial unentangled strictly asymmetric state |'i/'(0)) = 
|e,2)(i)|5,2)(2)|5,0)(3), that is, I?(0) = 1 in (C7), which 
lies in sector b. The state of the system is given for later 
times by the entangled state 

\i^{t)) = A{t) [|e, 0)(i)|<7, 2)(2) + |e, 2)(i)|g, O)^^)] |g, 2) (3) 



\(3) 



+ B{t) |e,0)^^^|g,4)^^) + |e,4)(^^|g,0)^^^ [5,0)^^^ (22) 



-C(t)|e, 0)(i)|g, 0)(2)|5, A)^^^+D{t)\e, 2)(i)|g, 2)(2)|g, 0)(3), 



where the probability amphtudes in (22) are given in Eq. 
(E3). Note that the system returns to its initial state 
for times r such that ^t = mr for n = 1,2,3,- •• and 
that D{t) y^ for t real. In addition, for any probability 
amplitude X{^t), \Xi^t)\^ = \X{t: - ^t)\^. 

Maximal entanglement is obtained by minimizing 
|C(t)p + |I?(t)P, which cannot vanish owing to the initial 
state being strictly asymmetric. Minimization occurs for 
^T w 0.1930,0.8542, where the minima have the same 
values and is |C(r)p + \D{t)\'^ w 0.1829. The measure 
of entanglement is log2(5.467) ~ 2.450. 

There are several simpler tripartite entangled states 
that occur at different times when some probability am- 
plitudes vanish: (a) ^Ti = 7r/5, B{Ti) = 0, C{Ti) = 0, 
|^(Ti)|2 = |sin2(7r/5) « 0.1536 and \D{Ti)\^ = 1 - 
|sin2(7r/5) « 0.6929 and so 2\A{Ti)\^ + \D{Ti)\^ = 1; 

(b) ^T2 = 7T/3, A{T2) = 0, B{T2) - 0, \C{T2)\'' = 18/25 
and \D{T2)\'^ = 7/25 with |C(T2)|2 + \D{T2)\'^ = 1; and 

(c) ^Ta = f - iarccos([V5 + l]/4) w 1.2566, B{T^) = 0, 
C(T3) = 0, |^(T3)|2 = |[1 + cos(7r/5)] « 0.4020 and 
\D{T2)\'^ = I - |cos(7r/5) w 0.1960 with 2|^(Ti)|2 + 

\D{T,r = 1. 
Is is interesting that for the above case with ^r = 7r/3, 

one has that A{t) = B{t) = with C(7r/3) = \/6(3 - 
V3i)/10 and D{T:/i) = (2 + V3i)/5. Therefore, the prob- 
ability that the system is in its initial state with two 
photons each in cavities 1 and 2 is 28% and that all four 
photons are in cavity 3 is 72%. 



y/l 1 



^^i:^b.8)'"te-»>'"i».») 



(3) 



(24) 



Note that if K{t) = 0, that is, ^r = 7r//\/66, / ^ 
0, 1,2,---, then F{t) = according to (D8); however, 
the converse does not follow. Therefore, for times ^r = 
{I + l/2)7r/V66, I = 0, 1, 2, • • •, F{t) = but K{t) = 
-2V3U/II, and so 



IV'(t)) 



(2)1 



-|5,2)W|5,2)(^)|g,2) 



v(3) 



2V30 1 



lf7!5:^l'?'6)^'^l5'0)^^^l.9,o) 



(3) 



(25) 



The geometric measure of entanglement of the state 
(25) is E{\iP{t))) = log2l21 « 6.92. This type of 
entanglement encoding is that of a W-state, a tripar- 
tite state of three qubits [33|. On the other hand, if 
T = {1 + l/2)7r/2%/66e, Z = 0, 1, 2, • • •, then 

|^(r)) = A|5,2)W|5,2)(2)|5,2)(3) 



^(-l)'7iE^l5^4)W|5,2)(^)|5,0) 



(3) 



(26) 



C. Initial state Eq. (21) 

In Appendix D, the explicit solution is given for an 
initially unentangled, symmetric state in all three cavi- 
ties, which can occur only in the A^ = 6 manifold, viz., 
the initial state given by Eq. (21) with a = c = e and 
b = d = f. Solutions (D8)-(D11) are not periodic owing 
to the incommensurate nature of the eigenfrequencies. 
However, for the initial state \g,2)(^^\g,2)^^^\g,2)'^^\ the 
solution is given by (D8), which are periodic with period 
T = 7r/\/66^ and so 

|V(t))=A(i)|.9,2)W|5,2)(2)|g,2)(3) 
+ ^W7fE^l5'4)*'^l5,2)(^)|5,0)(^) (23) 



X(t)i=^P|5,6)W|g,0)(2)|g,0) 



(3) 



V3 



where the probability amplitudes A{t), F{t), and K{t) 
are given by (D8) with a — 1. 

The amplitude A{t) vanishes for cos(2-\/66^T) — —5/6 
resulting in the maximally entangled state 



|^(r))=±^i=^F|g,4)W|5,2)(2)|.9,0) 



VgVg 



(3) 



^i=^P|g,6)W|g,0)(^)|5,0) 



(3) 



The geometric measure of entanglement of the state (26) 
121 
2^5" 



is E{\^{t)}) = log, i# « 2.28. 



VI. SUMMARY AND CONCLUSIONS 

We have studied the dynamical behavior of a deter- 
ministic system constituted by three identical cavities 
each enclosing a three-level atom with intracavity inter- 
actions governed by two-coherent-photon hopping. We 
consider atom/photon states for each cavity that cor- 
responds to a multileveled system of three-, five-, and 
seven-dimensional spaces. However, the dynamical state 
of the system, owing to conservation laws, lie in sub- 
spaces given by 6, 18, and 38 dimensions rather than the 
usual 33 = 27, 53 = 125, and 7^ = 343, respectively. 
Explicit analytic solutions are found for arbitrary initial 
unentangled pure states in the manifolds N—2, 4, and 6 
when the exchange of photons between cavities occurs at 
a much faster rate than the rate of atomic transitions, the 
large hopping limit. Tripartite entanglement between the 
cavities can be used to generate bipartite entanglement 
between any two cavities via local quantum operations, 
which can be performed either on the atom or on the 
photon pairs or even on both the atom/photon state of a 



given cavity. Dynamically generated tripartite entangled 
states, which reflect the symmetry of the Hamiltonian 
and the symmetry properties of the initially unentangled 
state, are given by the superposition of symmetric, anti- 
symmetric, and asymmetric states. The symmetric states 
are generalizations of the W-states. 



Appendix A: N=2 manifold 

The Hilbert space for three cavities with the state of 
each cavity given by the qutrit \g,0), |g, 2), and |e,0), 
normally lies in a Hilbert space of 3'^ = 27 dimensions. 
However, the dynamics of the N = 2 manifold is governed 
by only a 6-dimensional subspace owing to the constancy 
of the operator N with general state vector 

|V(t))-A(i)|.9,0)«|5,0)(2)|g,2)(3) 



+B{t)\g, 0)(i)|5, 2)(2)|g, 0)(3)+C(i)|5, ^)^'^\9, 0)^'^\9, 0)^''> 

(Al) 



(i)|p n\(2)u n\(3) 



+D{t) \g, Oy'> \g, 0) ('>\e, 0)^'>+E{t) \g, 0) ^^^ |e, 0)^^^ \g, 0) 

+F(i)|e,0)«|g,0)(2)|5,0)(3). 
The dynamical equations of motion are 

iA = A + 2^B + 2^C + tan 0oD, 



iB = 2^A + B + 2^C + tSin6oE, 



iC = 2^ A + 2^B + C + tan 9oF, 



D = tan 6*0 A + tan^ 6*0-0, 



E = tan OoB + tan^ OqE, 



(A2) 



A{0) + 2B{0)-C{0) 



^^2^(.t 



(A3) 



C{t) « - [A(0) + i?(0) + C(0) 



-Hit 



1 r 
3 



A(0)-B(0) + 2C(0) 



„2iCt 



D{t) w D{{)) E{t) w E{id) F{t) « i^(0), 

where the amplitudes are periodic with period T = tt/^. 
Results (9) in the text follow from (Al) and (A3) for 
A(0) == a, B(0) = 0, C(0) = 0, ^(0) = 6, E{0) = 0, 
and F(0) = 0. Note the difference in the labeling of the 
probability amplitudes in (Al) and Eq. (6). The average 
time the photons spend in each cavity depends on the 
initial state of the system and is bounded for cavity 3 
by i/J" \A{t)\^ ^ I + 5l^(0)P with analogous bounds 
for the amplitudes B{t) and C(t), for cavities 2 and 1, 
respectively. The equality holds when the photons are 
initially in cavity 3, that is, A{0) = 1, in which case the 
photons spend 5/9 of the time in cavity 3 and 2/9 each 
in cavities 1 and 2. 



Appendix B: N^4 manifold 

The Hilbert space for three cavities with the state of 
each cavity given by the five-dimensional qudit |5,4), 
\g,2), |5,0), |e, 0), and |e, 2), normally lies in a Hilbert 
space of 5'^ ~ 125 dimensions. However, the dynam- 
ics of the A^ = 4 manifold is governed by only an 18- 
dimensional subspace with general state vector 

|^(t))=A(i)|5,4)W|5,0)(2)|5,0)(3) 



+ l3,2) 



(1) 



B(i)|5,0)(2)|5,2)(3)+C(t)|.9,2)(2)|5,0)(3) 



where we have introduced the dimensionless time [{E^ — 
Eq) cos^ 6'o]t/7i -^ t and the dimensionless hopping cou- 
pling H/[{E^ - Eo)cos^ Oo] ^ ?. 

In the limit of large hopping strength ^ 3> 1, the solu- 
tions are 



A{t)^^\A{0)+B{0) + C{0) 



-Hit 



+i?(t)|e,0)(2)|5,0)(3)+iJ(t)|g,0)(2)|e,0)(3) 



1.9,0) 



(1) 



F{t)\g,0)^'>\gAr>+G{t)\gAr>\g,0) 



(2)u n\(3) 



+J(i)|e,2)(2)|g,0)(-^)+i^(t)|g,0)(2)|e,2)(3) 



2A{0) - B{0) - C(0) 



,2*5* 



B{t) 



A(0) + B(0) + C(0) 



-Hit 



L(i)|e,0)(2)|e,0)(3) +M(i)|e,0)(2)|5,2)(3) (Bl) 



-A^(t)|.9,2)(2)|e,0)(-^)+P(t)|.g,2)(2)|g,2)(3) 



|e,0) 



(1) 



iZ(t)|5,0)(2)|g,2)(3) + 5(t)|5,2)(2)|g,0)(3) 



-r(i)|e,0)(2)|g,0)(3)+C/(t)|5,0)(2)|e,0)(3) 



+W-(t)|e,2)(i)|g,0)(2)|5,0)(3). 

The dynamical equations for the probabihties aniph- 
tudes in (Bl), in the hmit of large hopping strength 
h^ ^ [{E2 — £'^)cos^02]: are given by the following 
uncoupled sets: For the amplitudes A, B, C, F, G and, 
P 



and 



lA « V24e(B + C), 

i{B + (7) « 2V2AS,A + 2S,{B + C) 

+ \/24e(G + F) + 4eP, 

i(G + /^) « V24e(B + C) + 2V24eP, 

iP « 2^(B + C) + V24e(G + F), 

i(S - G) « -2i{B - G) - y24^(G - F). 



(B2) 



(B3) 



i(G-F)«-V24e(F-G). 

the eigenfrequencies for the system of equations (B2) and 
(B3) are -8^, -6^, 4^, and 12^. 

For the amplitudes D, J, and M and 

it) K. 2£,J + 2CM, 

ij « 2^1) + 2£,M, (B4) 

iM w 2CD + 2iJ. 
For the amplitudes E, K, and A^ 

iF « 2^K + 2(N, 

ik w 2^E + 2^N, (B5) 

iiV « 2CF + 2^f:. 
For the amplitudes L, T, and [/ 

i « T « f7 « 0. (B6) 

Finally, for the amplitudes R, S, and VF 

iR « 2^5* + 2^H^, 

iS « 2Ci? + 2CH/, (B7) 

iW « 2Ci? + 2^5'. 

Solutions of the system of equations given by (B4), (B5), 
and (B7) are the same as solution (A3) for the system of 
equations (A2) in the limit ^ ^ 1. The eigenfrequencies 
for the Eqs. (B4), (B5), and (B7) are -2^, and 4^. 



Appendix C: N=6 manifold 

The Hilbert space for three cavities with the state of 
each cavity given by the seven-dimensional qudit \g,6), 
|5,4), \g,2), |g,0), |e,0), |e,2), and |e,4), normally lies in 
a Hilbert space of 7"^ — 343 dimensions. However, in the 
limit of large hopping strength h^ ^ [{E4 — E^) cos^ 64] 
the dynamics of the N = 6 manifold is governed by only 
a 38-dimensional subspace that in the limit of large hop- 
ping strength separates into four sectors: (a) all three 
cavities are in the ground state, (b) one cavity is in the 
excited state, (c) two cavities are in the excited state, 
and (d) all three cavities are in the excited state. 



1. Sector a 

The general state when all three cavities are in the 
ground state is 

|^(t))-A(i)|.9,6)W|.9,0)(2)|5,0)(3) 



+F(t)|g,4)(i)|g,2)(2)|5,0)(3)+G(t)|.g,4)(i 



+i?(t)|g,2)(i)|5,4)(2)|5,0)(3)-KF(i)|5,2)(i 



+F(t)|ff,2)(i)|5,2)(2)|5,2)(3)+G(t)|ff,0)(i 



-J(t)|5,0)W|g,4)(2)|g,2)(3)+i^(t)|g,0)(i 



-F(i)|5,0)W|g,0)(2)|g,6)(3 



l3,0)(^)|g,2)(3) 

|5,0)(^)|5,4)(3) 
(CI) 

lff,6)(^)|5,0)(3) 
lff,2)(^)|g,4)(^) 



The equations governing the time behavior of the prob- 
ability amplitudes break up into two groups, 

i(B - G) « 2^{C -B) + 12^(D - E) 



i{D - F) w 12^(F - G) + V60^(G - L) + 2^{J - K) 

i{G -L)k ^MS,{D -E + J-K) (C2) 

i{j -k)K 2S,{D - F) + V60e(G - L) + 12^{K - J), 
with eigenfrequencies 14^, —2^, and (1 ± \/241)^; and 
i{B + C) w 2V60eA + 2^(F + G) + 12C(F + F) + 2\/24^F 



i{D+E) « 12e(F+G)+V60^(G+F)+2^(J+F)+2V24CF 



z(G + F) « \/60e(F + F + J + F) (C3) 
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i{J+K) w 2^{D+E)+V60C{G+L)+12^{J+K)+2V2^F 



2. Sector b 



iF « V24^iB + C + D + E+J + K) 



iA w V60^{B + C), 

with eigenfrequencies 0, 2^, (—1 ± -\/241)^, and (7 ± 
VolS)^. Note that the eigenfrequencies are incommensu- 
rable and so the system never returns to its initial state 
albeit it can get arbitrarily close to it. 

Consider the case with initial condition A{0) = 1 in 
(CI). Eq. (C2) implies that B{t) = C(i), D{t) = E{t), 
G{t) = L{t), and J{t) = K{t) and so 

|V(t))=A(i)|5,6)«|5,0)(2)|g,0)(3) 



-S(t)l5,4)(i) 



E{t)\9,2)^'^ 



|g,2)(2)|g,0)(3) + |5,0)(2)|g,2) 



(2)1. n^(3)-^u n\(2)i. zi\(3) 



|.g,4)(2)|^^^0^(.)^|^^0^(.,|^^^4^ 



(C4) 



-GWl5,0)« 



|5,6)(2)|5,0)(3) + |g,0)(2)|5,6)(3) 



The general state for one of the cavities to be in the 
excited state is 

|V'(i))-A(t)|e,0)W|5,2)(2)|5,2)(3) 



+B(t)|e,0)(i)|5,4)(2)|g,0)(3)+C(t)|e,0)(i)|g,0)(2)|5,4)(3) 

(C7) 



(i)u n\(2)u 9\(3) 



+i?(t)|e,2)(^^|g,2)^2,|^^0^(.)^^(^)|g^2)^^)|g,0)^^)|g,2) 



+ni)|e,4)W|g,0)(2)|5,0)(3). 

The equations governing the time behavior of the prob- 
ability amplitudes break up into two groups, 

i{A -b)^ 2^{D -A) + ^/2A^{C - F) 



i{C-F)^V24CiA-D), (C8) 

with eigenfrequencies —6^ and 4^, and 
i{A + b)Ki 2^{A + D)+ 4CE + V24e(C + F) + 2V24^B 

i{C + F) w V24C(A + D) + 2V2A^E (C9) 



-if(i)|.g,0)« 



|5,4)(2)|g,2)(3) + |5,2)(2)|5,4)(3) 



+F(i)|.g,2)«|.g,2)(2)|.g,2)(3), 

with normalization |^(f)|2 + 2(|B(f)|2 + |S(t)|2 + G'(f)|2 + 
|i^(t)H + |F(t)|2 = l 

The system of equations (C3) can be solved explicitly, 
in particular, the probability amplitudes A(t) and F(t) 
for the unentangled states in (C4) are given by 



A{t) 



2 10 

\ ( 

11 29 



-2i^t 



66 



(1 



7 



/313' 



^{-7+V3T3Mt 



iB « V24^(v4 + D) 



lE « 2(,{A + D) + V24^(C + F), 
with eigenfrequencies —8^, —6^, 4^, and 12^. 

3. Sector c 

The general state for two cavities to be in the excited 
state is 

|^(i))=A(i)|e,2)«|e,0)(2)|g,0)(3) 



+ 7^(1 



7 



66^ ^313^ 



^-(T+ysTs)*^*^^!^! 



14, 8 



87^ 7V24T^ 



Jl+V24T)j5t 



(C5) 



14 



(1 



,(l-V24T)j{t 



and 



Fit) = 



87' 7V24T' 



10 VTo 



11 22 



1- 



/313 



-(7-y3T3)i{t 



22 



/3T3 



-(7+7313)44* 



(C6) 



+i3(i)|e,0)W|e,2)(2)|g,0)(3)+C(i)|e,0)«|e,0)(2)|.g,2)(3). 

(CIO) 
The equations governing the time behavior of the prob- 
ability amplitudes are 



iA « 2S,B + 2CC, 



iB w 2^A + 2^C, 



iC « 2^yl + 2e,B. 



(Cll) 



Solutions of the system of equations given by (C8) are 
the same as solution (A3) for the system of equations 

(A2). 
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4. Sector d 

The general state for three cavities to be in the excited 
state is 

|V(t))=A(t)|e,0)W|e,0)(2)|e,0)(3). (C12) 

The equation governing the time development is given by 

i«0 (C13) 

and so A{t) « A{0). 

Appendix D: Symmetric N=6 

One can obtain the 38-dimensional vector space for 
the TV = 6 manifold from (Al) by applying the following 
replacements: |(/,4)'^^^ by \g,6)^^\ \g,2)^^^ by a hnear 
combination of the vectors |g,4)(^) and \e,2)^^\ \9,0)^^'' 
by a linear combination of the vectors \g,2)^^' and 
|e,0)'-^'', |e, O)*^^' by a linear combination of the vectors 
|5,4)(i) and |e,2)(i), and |e,2)(i) by |e,4)(i). The 
preceding gives rise to 26 vectors. In addition, one must 
add to that set the 10 vectors that are the product of the 
vector \g,Oy^> and the linear superposition of the vec- 
tors |g,6)(2)|g,0)(3), |g,0)(2)|5,6)(3), I5, 4)(2)|e, 0)(3), 
|e,0)(2)|5,4)(3), |5,4)(2)|5,2)(3), |.g, 2)(2)|.g, 4>(3), 

|e,2)(2)|5,2)(3), |5,2)(2)|e,2)(3), \e,A)(^^g,0)(^\ and 
|g, 0)^^''|e,4)'^'^'. Finally, we have to add unentangled 
symmetric vectors, which cannot occur in the N — 2 
and TV = 4 manifolds, viz., I5, 2)(i)|5,2)(2)|c,, 2)(3) and 
|e,0)(i)|e,0)(2)|e,0)(3), for a total of 38 vectors. 

In the A^ = 6 manifold, the initial unentangled state 
that is totally symmetric in all three cavities is given by 

|V(0)) = \a\g, 2)(i) + 6|e, 0)^^^] \a\g, 2) (2) + b\e, O)^^) 



a|g,2)(3)+6|e,0)(3) 



(Dl) 



where 



|6p = 1. The case with 6 = reduces to 



that considered in Appendix C in sector a for ^ 3> 1. 
The trivial case a — gives a constant amplitude as 
indicated in (C13). 

The general symmetric state is given by 

|^(t))=^(i)|3,2)«|.g,2)(2)|g,2)(3) 






+i?(i)|e,0)(i)|e,0)(2)|e,0)(3) 



(3) 



(3) 



^W7fE^l3'4)(')|e,0)(2)|.g,0)(3) (D2) 



+^W7gE^l5'4>'''l.9,2)(2)|g,0) 



(3) 



"GW^E^I^'2)W|.g,2)(2)|g,0)(3) 



-^W7|E^I^'2)''^|e.0)(2)|g,0) 



(3) 



-K{t)^J2P\9,Q)^'^\9,0) 



(2) 



5,0) 



(3) 



J(t)i=5:P|e,4)W|g,0)(^)|5,0)(3), 



V3 



where for the amplitudes B{t), C{t), K{t), and J{t) 
the sum is only over even permutation of the cavities; 
whereas, for E{t), F{t), G{t), and H{t) the sum is over 
both even and odd permutation. The state (D2) brings 
forth all 38 states that span the A^ = 6 manifold, which 
corresponds to 10 symmetric, 4 antisymmetric, and 24 
asymmetric states that treats all the three cavities on the 
same footing. The antisymmetric states follow from the 
states associated with the amplitudes E{t), F{t), G{t), 
and H{t) in (D2), where a —1 is inserted before those 
terms arising from an odd permutation. 

The symmetric states are either entangled or unen- 
tangled, with the former coming in two varieties, viz., 
involving three or six states. The presence of only two 
unentangled symmetric states indicates that in the A = 6 
manifold one does not have the general GHZ maximally 
entangled state for the seven-dimensional qudit associ- 
ated with each of the three cavities in the space spanned 
by 7"^ = 343 vectors. 



|Gi?Z> = i=^b-)Wb-) 



(2)|^y3) 



(D3) 



where the sum over j is over the states \g, 6), \g, 4), \g, 2), 
\g, 0), |e,0), |e, 2), and |e,4) and the trace over one of the 
three cavities gives rise to an unentangled mixed state. 
In the A^ = 6 manifold, the trace over one of the three 
cavities of the symmetric states involving six states pro- 
duces an entangled reduced density matrix for each cav- 
ity given by the three states that span that space. How- 
ever, a trace over two of the three cavities, for instance, 
the reduced density matrix for each cavity associated 
with the probability amplitude E{t) in (D2) is given by 
p=i[|g,0)(g,0| + je,0)(e,0| + |5,4)(5,4|. 
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In the large hopping Umit ^ ^ [{E^ — E^) cos^ 64] the 
equations for the probabihty amphtudes decouple and 
are as follows: For A{t), F{t) and K{t) 

dA 



.^«12^A- 



2\/30^X, 



(D4) 



at 



for B(t), £;(t), G{t), and J(t) 



dt 



dE 
'^~dt 



4V3^E + 2V2^G, 



iVS^B + 2V6^G, 



(D5) 



dG 
i Si 2V2CB + 2V6^E + 4\/3C J, 
dt 



for i:>(t), 



and for C{t) and iJ(i) 



.dJ 

l—- K, 

dt 


4V3eG, 


.dD 

J— — 
dt 


-«0, 


dC 

i—— PS 
dt 


2V2^H, 


dH 


2V2£_C. 



(D6) 



(D7) 



The initial conditions for the state (Dl) correspond 
to A(0) = a^, B{0) ^ VSa^b, C(0) = VSab^, and 
D(0) = 6'^. Therefore, the probability amplitudes in 
the different sectors are separately conserved and so 
\Ait)\' + |F(t)|2 + \K{t)\' = \a\^ \Bitr + \E{t)\^ + 
|G(t)|2 + |J(i)|2 = 3|an5|2, \D{t)\^ = |6|6, and \C{t)\^ + 
\H{t)\'^ = 3|ap|6|^, which when summed gives the over- 
all probability of unity. Note that if a = 0, then 
the three cavity system remains in its initial state viz., 
|e,0)(i)|e,0)(2)|e,0)(3). However, if 6 = 0, viz., the ini- 
tial state is \g, 2)^^'>\g, 2)^'^^g, 2)^^\ then the only nonzero 
amplitudes are A{t), F(t), and K{t) with normalization 
condition |A(f)|2 + |F(t)|2 + \K{t)\^ ^ 1. 

The general solution of the system of equations (D4)- 
(D7), with initial condition (Dl), are for A{t), F{t), and 
K{t), 



A{t) 
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6cos(2V66^i) 



F{t) = -—a^V66ism{2V66&), (D8) 



K{t) = — VaOa'^ 



cos(2V66^i) - 1 



For B{t), E{t), G{t), and J(i), one has exact solutions; 
however, we present their simpler numeric rather than 
the purely analytic results and so 

B{t) = a26[0.4054e-"-2644««* + 0.3995e-3-^306^?* 



-H0.0838e' 



,8.6745iCt 



0.8433e' 



6.3205iCt 



E{t) = a% 0.4607e-"-2644^«* + 0.3401e-3-^306»«t ^g) 



-0.2040e«-6745^et _ o.5968e6-3205.«* , 



G{t) = a^b 



0.4860e-"-2644*«* - .3061e-3-^306»ct 



-h0.2427e 



8.6745iCt 



.4227e' 



6.3205iCt 



J{t) ^ a% 



0.2989e-"-2644^«* - 0.5684e-3.7306«ct 



-0.1939e«-6^45«ct ^ o.4633e6-3205^?* 

For D{t) one has that, 

D(t) = b^. 
Finally, for C{t) and H{t) one has that, 
C{t) = ^/3a&2cos(2^/2^i), 

i7(t) = -\/3a62isin(2V2^t). 

Appendix E: Asymmetric N=6 



(DIO) 



(Dll) 



In Appendix D, the case of a symmetric initial state in 
all three cavities (Dl) was considered in the limit of large 
hopping strength ^ 3> 1. In this Appendix, we consider, 
in the same large hopping limit, the case where the atom 
is in the initial unentangled asymmetric state |V'(0)) = 
|e,2)(i)|g,2)(2)|c,,0)(3), viz. £1(0) = 1 (see (El) below). 
Note that in the limit of large hopping, only photons are 
transferred between cavities with no "transfer" of atomic 
excitations; therefore, the atom in cavity 1 remains in 
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the excited state and so the state of the system at later 
times is given by 



(l)|n 9\(2)U 9\(3) 



|V(t)>=A(<)|e,0)W|5,2)(^^|.9,2) 



+ |e,0)« 5(i)|.9,4)(2)|g,0)(3)+C(t)|g,0)(2)|.g,4) 



+ |e, 2)(i) [D{t)\g, 2}('^\g, O)^^) + E{t)\g, 0)^%, 2) (3) 

(El) 

+F(i)|e,4)«|g,0)(2)|5,0)(3). 

The equations that govern the time dependence of the 
probabihty ampUtudes in (El) are 

dA 

i—r- 
dt 



24^5 + V24CC + 2^D + 2^E, 



dP 
i r^ y/24(,D + V2A^E. 
dt 

The state (El) can be written as the sum of symmetric 
and antisymmetric states under the exchange 2 <-^ 3. The 
probabihty amphtudes A{t), B{t)+C{t), D{t)+E{t), and 
F{t) are the symmetric amphtudes with eigenfrequencics 
-8^, -6^, 4^,' and 12^. The amphtudes B{t) - C{t) and 
D{t) — E(t) are the antisymmetric amphtudes with eigen- 
frequencics —6^ and 4^. Accordingly, (E2) can be written 
as a set of two and four uncoupled equations for the an- 
tisymmetric and symmetric amplitudes, respectively. 

The solution of (E2), for the initially uncntangled state 
|V(0)) - |e,2)W|g,2)(2)|p,0)(3), that 'is, D{0) - 1, is 



A{t)^E(t)^^ 



. 2e-4««* - 3e6*«* + 2e*^*«* + 3e-^^'^' 



B{t) = F{t) 
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-4i4t _ „6iCt _ pSJ^t 



-12i{t 



(E3) 



i « V2i^A + V24^i:», 

dt 



m-^ 



2e-4^«* + 2e6^«*-e*^*«*+e-i2»«* 



dC 
dt 



(E2) 



D{t) 
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4e-4^«* + 6e6^«* + 2e**«* + ie-^^'^' 



i ~ 2CA + \/24eS + 2S,E + \/2Ae,F, 
dt 



i—-K2S,A + \/24^C + 2^D + \/24^F, 
dt 
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